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The ultrahigh energy neutrino cross section is well understood in the standard model for neutrino 



energies up to 10 GeV. Tests of neutrino oscillations [y^ 



from extragalactic sources of neu- 



trinos are possible with large underground detectors. Measurements of horizontal air shower event 
rates at neutrino energies above 10 10 GeV will be able to constrain nonstandard model contributions 
to the neutrino-nucleon cross section, e.g., from mini-black hole production. 



I. INTRODUCTION 



of neutrinos between the three neutrino flavors 1 
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The existence of ultrahigh energy (UHE) cosmic rays 
comprised of nucleons guarantees the presence of ultra- 
high energy neutrinos. Over megaparsec distances, UHE 
nucleons interact with the 2.7K microwave background 
radiation to produce the delta resonance, which decays 
to a pion and a nucleon. Charged pions decay to neu- 
trinos and muons, followed by muon decay, leading to 
fluxes of electron and muon neutrinos. Beginning with 
the pioneering work of Ref. 1, and continuing to recent 
improvements^ 01 j there are quantitative predictions for 
fluxes of neutrinos with E v ~ 10 8 GeV and higher ener- 
gies. 

At somewhat lower energies, one expects that sources 
of TeV photons such as active galactic nuclei should also 
produce neutrinos. The attenuation length of the pho- 
ton effectively cuts off the photon energy at Ey ~ 1 TeV 
for sources farther than 100 Mpc from the Earth. The 
long interaction length of neutrinos means that neutrino 
attenuation from distant sources is not an issue over as- 
tronomical distances due to low densities. 

A common feature to neutrino flux models is that 
roughly two muon ncutrinos/antincutrinos appear for 
each electron neutrino/antineutrino at the production 
site. (In the following discussion, we do not distin- 
guish between v and v.) In the past few years, the Su- 
perKamiokande experiment 3 together with the SNOQ 
results point to physics beyond the standard model: neu- 
trino masses and mixing. Parameters in the mixing ma- 
trix lead fairly robustly to an equal distribution of fluxes 



F = F = F 



(1) 
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after propagating over astronomically large distances. 
UHE tau neutrino interactions with nucleons in the Earth 
and atmosphere may provide signals of this phenomenon. 
Tests for deviations from the standard model ultrahigh 
energy cross section are another probe of non-standard 
model physics with UHE neutrinos. 

The neutrino energy regime determines the detection 
method. For E v < 10 6 GeV, underground detectors 
can be used to record upward muons and electromag- 
netic/hadronic showers. Experiments include water or 
ice Cherenkov detectors like AM AND A 0, NESTOR[|, 
ANTARES0 and IceCube^. Attenuation of the neu- 
trino flux in Earth is important, especially the difference 
between and v T attenuation|lll 1 1 2L Il3| . In the en- 
ergy range above 10 8 — 10 9 GeV, detection of horizon- 
tal air showers and radio detection are used to deter- 
mine the UHE v fluxes and to probe the neutrino cross 
section [Til 

13 El 13 13 E3 At intermediate energies, 
the double-bang signal of both production and decay ver- 
tices of the v T — ^ r — » v T X process can be used to iden- 
tify i/ r s|3|. Other detection possibilities include using 
the Earth as a tau neutrino converter, then searching for 
the shower from the decay of the emerging tau[21|. 

We present results in two of these examples: upward 
neutrino fluxes and their detection, and the potential for 
limits on the non-standard neutrino cross section from 
the OWL[l6j proposed detector. In this latter case, the 
non-standard model considered is mini-black hole pro- 
duction and decay [3 El EE 13, E3, a fea- 
ture of models with TeV-scale compactification of extra 
dimensions 29| . Before going into detail of these two ex- 
amples, we review the theoretical status of the neutrino 
cross section at high energies. 
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II. UHE NEUTRINO CROSS SECTION 

The neutrino-nucleon cross section is straightforward 
to calculate in the standard model 30]. The charged cur- 
rent cross section is 



q(x,Q) + (l-y) 2 q{x,Q) 



dx dQ 2 




(2) 

for W boson momentum squared q 2 = —Q 2 , lepton en- 
ergy difference g° = E — E' = u, Bjorken x = Q 2 /(2Mv) 
and y — v/E with nucleon mass M and W mass M\y- 
The parton distribution function combinations q(x, Q) 
and q(x,Q) are detailed in, for example, Ref. 31. The 
cross section depends crucially on the parton distribution 
functions |3^. As m(Q 2 ) gets large, the parton distribu- 
tions at small- x increase, however, the W boson propa- 
gator decreases with increasing Q 2 . At sufficiently high 
energies, the value of Q 2 saturates at Q 2 ~ M^, so the 
relevant values of x for a given incident neutrino energy 
E is approximately 



Ml 

2ME 



(3) 



This translates toi~ 1CT 2 for E ~ 10 6 GeV, and x ~ 
1(T 8 for E ~ 10 12 GeV. 

Measurements of the parton distribution functions 
have been made down to x ~ 10~ 6 at Q 2 ~ 0.1 GeV 2 
at HERA[33[. At a more relevant Q 2 ~ M 2 F , parton 
distributions have been measured to x ~ 10 _3 [34[. 

Beyond the measured regime, one is guided by theo- 
retical arguments concerning the small- a; parton distri- 
bution functions and their QCD evolution. As discussed 
in detail in Ref. 35, for xg(x, Qo) ~ x~ x at x <C 1, then 
xg(x, Q) ~ x~ x for A > 0.3. At small values of x, the 
quark distribution functions are dominated by the quark 
sea, which come from g — > qq, so the sea distributions 
have the same power law dependence: xq(x, Q) ~ x~ x . 

Modern parton distribution function fits have A ~ 
0.3 — 0.4, and extrapolations to small- x by this power law 
agree well with perturbative QCD evolution, as demon- 
strated, for example, by Gliick, Kretzer and Reya in Ref. 
36. Kwiecinski, Martin and Stasto have used a BFKL- 
type evolution at small- X and obtain cross sections that 
are in agreement with power law extrapolations |3?j . 

At very small values of x, recombination processes such 
as gg — > g should become important, ultimately leading 
to a flattening of the neutrino cross section as a function 
of energy. The recombination regime is unimportant for 
neutrino energies up to 10 12 GeV, though it should be- 
come important at even higher energies |38|. 

Fig. 1 shows the neutrino interaction length, propor- 
tional to the inverse of the cross section: 
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(4) 



with Na equal to Avogadro's number. For neutrino en- 
ergies above a few tens of TeV, the Earth's diameter is 



Earth diameter 



"Moon diameter 



I Horizontal atmospheric depth: 
'° r 3.6 x 10 4 cm.w.e 
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FIG. 1: The standard model neutrino interaction length in 
cm.w.e. = g/cm 2 , with reference distances through the Earth 
and the Moon. 



larger than the attenuation length. Even for neutrinos 
with E = 10 12 GeV, the horizontal atmospheric depth is 
short compared to the neutrino interaction length. 



III. UPWARD NEUTRINOS IN 
UNDERGROUND DETECTORS 

In our discussion of upward neutrinos, we restrict our 
attention to neutrino energies below 10 6 GeV 12]. For 
the upward neutrino flux, depending on the nadir angle, 
attenuation due to neutrino passage through the Earth is 
important. Neutrino oscillations of <-* v r lead to the 
arrival of equal fluxes of f M and v T at the surface of the 
Earth, however, it has been noted that the attenuation of 
the two species is quite different due to the short lifetime 

of the t|ijLIi1Ii3- 

Muon and tau neutrinos have the same neutral current 
and charged current cross sections above 1 TeV. (At 1 
TeV, the charged current cross section for v T — > t is 
about 5% lower than the charged current (CC) cross 
section due to the r mass. At 100 GeV, the v T CC cross 
section is about 25% smaller 39].) The difference in the 
attenuation of and v T fluxes has to do with the fact 
that muons lose energy as they traverse a medium before 
they decay. Electromagnetic energy loss for muons is 
significant, so the neutrinos from muon decay have much 
lower energies than the parent muon energy. On the other 
hand, electromagnetic energy loss for r's is important 
only for tau energies above ~ 10 8 GeV^(J. The muon 
neutrino and tau neutrino transport equations are 



dF^(E.X) _ F„JE,X) 



dX 



(5) 



dE y G u ^ 



u *(E,E y ,X) 
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FIG. 2: The ratio of the attenuated flux to incident flux for 
v T (upper curve) and (lower curve) for F u ~ E~ 2 at nadir 
angle 0. 



FIG. 3: The ratio of the attenuated flux to incident flux for 
v T (upper curve) and (lower curve) for F° ~ _E _1 at nadir 
angle 0. 



dF„ T (E,X) _ F Vr (E,X) 
dX Cl nt (E) 

POO 

+ / dEy [ff^^ {E, E y , X) + G™ {E, Ey, X) } 
JE 

8F T (E,X) F T {E,X) 



dX 



+ / dE y G Vt ^ t (E, Ey,X) 
JE 

where, for example, 
G v ^{E,E y ,X) = 



F Vt (E v ,X) 



dn 



NC 



dE 



-(E y ,E) . 



in terms of the neutrino flux F Vt and the cross sec- 
tion normalized neutral current energy distribution 
dn NC /dE. We have made the approximation that tau 
energy loss is negligible, and that only the decay of the 
tau is relevant, neglecting tau CC interactions. 

The feed-down of neutrinos from v T — > r is evident 
in comparison of the attenuated v r flux with the 
flux[ll[l3|. In Figs. 2 and 3 we show the ratio of the 
attenuated flux to the incident flux for v T (upper curve) 
and (lower curve) at nadir angle for two different 
fluxes, parameterized by 



F° = £?-t-V(1 + E/E Y 



(6) 



for E = I0 8 GeV. In Fig. 2, 7 = I and a = 0, and 
in Fig. 3, 7 — and a = 2. The tau neutrino pile- 
up is more pronounced for the flatter fluxes. At larger 
nadir angles, the pileup reduces because the number of 
interactions reduces, but the attenuation is also less. 

Experimentally, one would like to exploit the effect of 
the tau neutrino pile-up and to identify the tau neu- 
trino component of the flux. In a manner similar to 
the lower energy comparison of the neutral current to 



charged current event rates for neutrinos with ener- 
gies above 1 TeV, one can compare the upward electro- 
magnetic/hadronic shower event rate from 



v T N -> 

V-r_„_ e N — > V 

e 



v P N - 



t + hadrons, r — > v T 
t + hadrons, r — > v T 
hadrons 



hadrons (7) 
• e + v P 



r./i,e 

hadrons 



with the upward muon event rate from 



fi - 

T - 



X 

x, 



(8) 
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We have made a detailed comparison in Ref. 12. By mea- 
suring the upward shower and muon event rates with an 
energy threshold of 10 TeV, one should be able to see the 
effect of the tau neutrino component of the incident neu- 
trino flux. This is relatively independent of the a priori 
unknown spectrum of the incident flux, except for the 
overall normalization (so that there are enough events). 

An interesting possibility for a class of input neutrino 
fluxes is the enhancement of the upward muon rate from 
v T — * t — * where the final muon neutrino interacts to 
produce a muon[4l|. The muon neutrinos from tau de- 
cays are called "secondary neutrinos." There is incremen- 
tal energy loss at each of the processes: CC production 
of the tau, then muonic decay of the tau. This results in 
a contribution to the muon neutrino flux 42] , as seen in 
Fig. 4 for F ~ 1/E for nadir angle 0. The upper dashed 
line is F v jF® for the v T flux, the lower dashed line is 
for excluding the secondary neutrinos, while the solid 
line is the secondary muon neutrino flux. The implica- 
tions for the upward muon rate are shown in Fig. 5. The 
secondary neutrino flux for the 1 / E 2 incident flux is rel- 
atively less important. For less steep fluxes than 1/E, 
one expects that the secondary flux will be enhanced. 
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FIG. 4: The ratio of the attenuated flux to incident flux for FIG. 6: The cosmogenic electron neutrino plus antineutrino 
v T (upper dashed curve) and (lower dashed curve) and flux of Engel, Seckel and Stanev. 3 



secondary (solid curve) for F° ~ E 1 at nadir angle 0. 
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FIG. 5: Upward muon event rate for F° ~ 1CT 13 (GeV/E) 
(cm 2 s sr GeV) -1 without secondary neutrinos (dashed curve) 
and with secondary neutrinos (solid curve) for E„ > 10 4 , 10 5 
GeV. 



IV. HORIZONTAL AIR SHOWER SIGNALS OF 
MINI-BLACK HOLES 

As Fig. 1 indicates, the standard model neutrino in- 
teraction length is large compared to the horizontal at- 
mospheric depth (zenith angle 90° at sea level). Un- 
like strongly interacting particles, neutrinos can pene- 
trate deep into the atmosphere. Horizontal air showers 
initiated deep in the atmosphere are a signal of neutri- 
nos. Non-observation of an enhancement of horizontal air 
shower event rates limit a combination of the neutrino- 
nucleon cross section and the neutrino flux |18lll9| . Using 
the cosmogenic neutrino flux from cosmic ray scattering 
with the 2.7K photons, it is possible to set limits on non- 
standard contributions to the cross section. 

The cosmogenic neutrino flux, for the sum of electron 



neutrino plus antineutrino, evaluated by Engel, Seckel 
and Stanev[2| (ESS) is shown in Fig. 6. Their results 
are shown for two models of source evolution: one pa- 
rameterized like (1 + z) 3 for redshift z < 1.9 (standard 
evolution), and the other scaling like (1 + z) 4 for the same 
redshift (strong evolution). 

Discussed here are contributions from mini-black holes, 
predicted by models with the Planck scale in the TeV 
range and a number of extra dimensions |2flj. and the po- 
tential of an orbiting observatory to detect the enhanced 
rate of horizontal air showers |43j. There is an exten- 
sive literature on black hole production and detection 
by neutrino telescopes[H El El El El El El El 
It appears that enhancements of the cross section from 
i-channel Kaluza-Klein graviton exchange is not an im- 
portant effect |H El- 

Black hole production in these models is determined 
by the geometrical cross section, depending on the 
Schwarzschild radius. The Schwarzschild radius rs, as- 
suming n compactified dimensions and Planck scale Mo, 
is given bv|i^ 



rs 



1 




(9) 



Here, Mbh is the mass of the black hole. The geometrical 
cross section for neutrino interaction with parton j is 
then 



a{vj BR) 



"l(M BH 



1U BH ) 



(10) 



for s equal to the neutrino-parton center of mass energy. 
Eq. (10) enforces the requirement that s be larger than 
some minimum energy-squared, certainly larger than the 
Planck scale. Convoluting with the parton distribution 
functions /j after setting s = x s give a non-standard 
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FIG. 7: The black hole production cross section for n = 4 
(solid) and n = 6 (dashed) for several values of the minimum 
black hole mass, for a Planck scale of Md = 2 TeV. The 
standard model cross section is also shown (SM) with the 
dot-dashed line. The limits of Tyler, Olinto and Sigl 18 (TOS) 
and Anchordoqui et al. 19 (AFGS) are shown with the dotted 
lines. 
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FIG. 8: The event rate per year for black hole production for 
2 monocular eyes for OWL from the ESS cosmogenic electron 
neutrino flux (standard evolution), for n = 2,4,6 extra di- 
mensions, M D = 1, 2, 3 TeV as a function of M^/M D . Also 
shown is the standard model event rate. 



model contribution to the neutrino-nucleon cross section: 

a(uN - BH) = L»»* )2 dx &f H {xa) h{x,Q 2 ) . 
i J — ^ — 

(11) 

The black hole is expected to rapidly evaporate into 
standard model particles, and we assume that all of 
the ener gy of the black hole goes into the shower it 
produces|44j. 

Experimental constraints on deviations from Newto- 
nian gravity exclude n = 1 and Md ~ 1 TeV 19]. For 
n = 2, Md > 3.5 TeV is consistent with tests of the 
inverse square law.|45| while for larger n, there are no 
constraints from tests of Newtonian gravity for Md ~ 1 
TeV. Collider experiments put a lower limit on Md at 
approximately 1 TeV [iH Erf. 

Fig. 7 shows a representative set of parameter choices 
for the number of extra dimensions n and the minimum 
black hole mass, given a Planck scale of Md = 2 TeV. 
Also shown are curves representing upper limits on the 
cross section from Tyler, Olinto and Sigl 181 (TOS) and 
Anchordoqui, Feng, Goldberg and Shapere[l!| (AFGS), 
were they have used cosmogenic fluxes similar to that of 
ESS. These upper limits come from the AGASA and Fly's 
Eye limits on an excess of horizontal air shower events. 

The largest enhancement of the cross section comes at 
the highest energies. Detectors optimized to this energy 
range will be able to put the most stringent limits on the 
cross section. The pro pos ed Orbiting Wide- Angle Light- 
Collectors Experiment Ha| (OWL) and the Extreme Uni- 
verse Space Observatory [l7j (EUSO) are two such detec- 
tors. The OWL experiment would involve two satellites 
with photodetectors, orbiting 640 km above the Earth. 
The EUSO experiment would be located on the Inter- 



national Space Station, 380 km above the Earth. The 
conversion factor from the 2 telescopes to one at a re- 
duced altitude amount to a reduction of the event rate 
by a factor of about 0.2 compared to OWL. 

The threshold for OWL detection of horizontal air 
showers is on the order of 10 10 GeV. By rescaling the 
effective aperture for OWL detection of electron neu- 
trinos by the ratio of the standard model plus black 
hole cross section to the standard model neutrino-nucleon 
cross section, one can arrive at event rate predictions 
given the ESS incident electron neutrino flux. The range 
of electron neutrino induced event rates goes from 100 
events per year for n = 6 and M D = Mgg 1 = 1 TeV 
to less than one event per year for Md = 3 TeV and 
Mgjj 1 = 10 • Md for the standard evolution of the cosmo- 
genic flux as shown in Fig. 8. The standard model event 
rate is about 0.5 events/year. The strong evolution of 
the cosmogenic flux gives a factor of 2 enhancement for 
both signal and background. The OWL rates are about 
a factor of 20 times larger than the IceCube downward 
contained event rates for an energy threshold of 10 s GeV. 
IceCube will optimally detect fluxes with large low energy 
components, e.g., F v - \/E 2 27\. 

If all three neutrino flavors are included in the event 
rate calculation, the standard model rates increase by a 
factor of about 2, while the black hole rates are increased 
by a factor of 3 since the black hole cross section is flavor- 
blind. (The standard model rate depends on how much of 
the incident neutrino energy is translated to the shower 
in the final state.) Including these factors, one finds that 
OWL would be able to probe the fundamental Planck 
scale up to 3 TeV for n > 4 even with one year of data 
taking. 
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V. FINAL REMARKS 

Described above are two ways to search for nonstan- 
dard model physics with ultrahigh energy neutrinos: 
searching for evidence of — > v T oscillations with tau 
appearance in upward muon and upward shower rates, 
and mini-black hole production in neutrino induced hor- 
izontal air showers. Nonstandard model physics may be 
manifest in the flux of neutrinos, through exotic particle 
decay to neutrinos |47|. A wide range of experiments: un- 
derground Cherenkov detectors, surface arrays and satel- 
lite light collectors cover a huge range of energies from 
GeV to 10 12 GeV incident neutrinos. The view from the 



ultrahigh energy neutrino window on astrophysics and 
particle physics should provide interesting tests of the 
standard model and its extensions. 
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